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A new transformation of the Klein integrands is 
introduced so as to avoid the unnecessary inconvenience 
of having the integration variable different from that pro- 
vided by experiment . As an alternative to the Rydberg-Uein- 
Rees equations, the simplification of Singh and Jrdn is ex- 
t endec to eliminate all reference to rotation and Alain’s 
troublesome g function not only do numerical integrations 
disappear, but *»‘K8 results are reduced to their bare essentials. 
Application of this extension to the fourteenth vibrational 
level of norral no I ecu la r hydrogen, which energy corresponds 
to 98 % of the way tcvnrd dissociation, reveals an error of 
only 0,41 % in the maximum turning point. 

IRYROt'tCTIC- : 
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The Ryiberg-Rloih-ilees technique, hereafter called RfR, is a method 
for constructing diatonic potential. envv-qy curves ~ d Ljfec.t ly free spectroscopic 
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tat* without the prior assumption of as analytic force law. ’Tenewed in- 

4 

terest i» the procedure occurred about 19S9 when Venders 1 iee , et al , 
applied the equations to hydrogen for the ultimate p^pose of calculating 
atom-atom collision cross sections sou high- temperature cas transport co- 
efficients; other computations followed i mediately and have been summarized 

by Zelexnik . S 

Perhaps the most distinguishing feature of much of this work, especially 
the earlier papers, is the persistent adherence to Sees'* integral evaluations 
for energy level s containing no higher than second-order vibrational contri- 
butions. Molecular states for which the experimental data could net be repre- 
sented over the entire range by a single quadratic were, instead, piecewise 
fitted over snail er intervals and correspondingly integrated. This procedure, 
being open to Question even at the lower levels, could yield quite substantial 
errors near the dissociation limit. The ot«Jy other strictly analytic approach 

of note is Jaraaia's scries development , ” wherein the set of vibrational levels 

7 

is expressed in the fora of an unrestricted polynomial and hunha&'s potential 
is inverted for the classical turning points. 

After a brief discussion of the validity of WKS and HP? methods, the re- 
maining purposes of this paper are three- fold as follows; to remove singu- 
larities from Klein integrands , to develop a substitute formulation in which 
such Integrals disappear from the working equations , an 5 to compare results with 
acre rigorous calculations in the case of ground- state molecular nydrogon near 
its dissociation limit. 
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l JAM POT'.aTIAL 


Nearly ail semi -empirical schemes for counting diatomic force laws 
originate from the slightly modified Bohr-Aewserfeld (first-order WK2) 


quanta® condition 



^(v+5) 
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where r is the equilibria* separation between nuclei and v is the vi- 
£ 

irrational Quantum number. It is possible to solve this equation either 
directly for the classical turning points of an arbitrarily chosen vibra- 
tional energy, as did Pyclberg and Klein, or by assuming various ssries ex- 
pansions for the potential anti its first derivative with respect to £. Using 


V * a Q i 2 (i+iij €+3, ^ + ) - C 2 ) 
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thmhiKs evaluated the phase integral and obtained 
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for the energy level G v . 

Since both the BX,& and Dunham expressions derive frm the WKh aunroxi- 
metion, and since that approximation regard's h as an expansion parameter, 
neither theory is a p rior i valid to first order except for high vibrational 
quantum numbers and correspondingly large displacements from equilibrium of 
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the classical turning points. It thus follows that sines the ccsyarat ive ly 
limited supply of such high energy levels say provide but a snail fraction 
of the number of coefficients in cq. (2) needed to precisely describe the 
true potential at th4£~ heights, the resulting truncated series may fit very 
poorly in ot£er~ Joeains of t. The danger would be especially great near 
the potential minima^ A >+*~ — *r (jr > 

GturjL+*ixct*.' 

Despite the coa&oa origin of Klein (J*KR) and Dunham potentials, scv«nl 

sis conceptions have developed sines the original work: It is often stated,** 

for example, that the semielcssi csl foundations of the P£F wethed render it 

untenable at low v, whorers >mha*j ! s first approximation is accurate near 

potential rdisima. Perhaps the difficulty arises from looking ton closely at 

the form of £q. (2) without recognizing the role played by WKB theory in 

determining th* coefficients as functions of spectroscopic constants. In any 

event, the distinction is somewhat academic for many aolec.Iar states because 

2 

corrections to (1) of order h and hipber are completely negligible; hence 
no con si deration of the c lassies 1 limit is necessary and all restrictions 
oa v isajr be removed. Accordingly, both, the Klein and Dunham potentials will 
suffice over the entire ;■ range in certain cases - but to the extent that either 
fails near the *dnirum, so also isust the other. 

X - : ' 

' $ A COlvYHNlEAT TRANSFORMATION 


The RKR equations Cot a nonfat at ing molecule are listed as fellows: 


2 .1 1/7 

t . * (f * i f) *f. 
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where G ana B are the pure vit rational an# first -oruer rotational 

v v 

energies, respectively, each measured in inverse centimeters, and h is 


Planck’s constant in erg-seconds. Other parameters include the reduced mss 
# . the speed of light c, and the vibrational puantus number v corres- 
ponding to that energy for which the classical turning noints T miiX 

r , are renal rad. 
win 

With the assumption now that G is an analytic function of v, the 
integrand of <7> easy be written 
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where the primes and superscripts on G all refer to differentiations of 
various orders with respect to v. Suggested immediately by ($) is the 

f /2 

simple transformation from v’ to x» (v-v * ) *' , according to which Klein’s 


f and 


g functions t*ecot« 
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with 


‘ *(V.x) » B - « (V-X 2 *ij + Y. (v-x‘ ♦ i) + 


( 12 ) 


Several advantages of e quation." (10) and (11) ever (7) and (3) are 


obvious: Me longer do integrand singularities exist, only the properties 

of the G y curve at the chosen level arc required, and x is a convenient 

variable of integration. This last advantage becomes evident when comparison 

10 1/2 

is made with ‘dels smart’s substitution ' of x for (C - G ) ' to obtain 

V V ■ 
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Numerical evaluation of (15) introduces an unnecessary inconvenience: Since 

the integrand is available in v' rather than x, only increments of the 
former may be selected arbitrarily and these of the latter then calculated 
from Ax - (G y - C yJ + ffl ) 1/2 - (G y - G y ,) i/2 . 

It is further noted that the upper limits to the summations in Bqs. (9)*“ 
( i ft, - 4 - 11) are, at most, the degree of the highest polynomial with which the 
spectroscopic data cn G y can bo fitted, i,e., the number of observed vi- 


brational levels. In addition, because of the relative unirnortance of 
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higher derivatives as compared with the first, sue.*; polynomials ordinarily 
■ay be considerably truncated by least squares methods and no significant 
error is expected froa the appearance of successive slopes of experimental 
measurements. 

Zelsodl^ has recently developed a nw method for evaluating the Klein 
integrals that alee essentially retains the original integration variable, 
he removed the singularities with an integration by parts and expressed the 
remaining integrands in terms of polynomial expansions containing derivatives 
with respect to v*. ;*© particular gain is apparent, however, over Eos. (10) 

and (11) 1 b which such derivatives at all points below the chosen v are 
specifically elisinated. Other, more approximate procedures involving ana- 
lytic curve fitting only near the upper limits of Hqs. (7) and (8) with 
numerical integrations over the rest have been applied by several investi- 
gators, as has the technique of Gaussian integration** near the singu- 

larities. 


THE MORSt POTE^mAL 


As a prelim nary to our alternative to the HK8 set hod , it is convenient 

lb 


to point out the simplifications inherent ia the Morse potential 
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the differentiation of which yields 
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for substitution into Dunham’s expression (5), Since 
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the solution of Eqn (5): may be written is closed for» as follows 


accordingly. 
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( 1 *) 


/V«c« j 7 

Davidson aiid Pines apparently wore the first to recognize that 
rigorour Morse energy levels arc identical (apart from half-integral quantum 
numbers) with those rrecicted by the 3ohr- '-eta^erf e 1 d quantum condition, a 
fact which also follows f^om the vanishing ox first and higher corrections 
to (3). In spite of the support it 'rents to funh.’ws and UKR methods, this 
agreement seers to hsve 1 ecu overlooked in the? marc recent literature, 

Verification of the above statement in the case of Klein's f function 
proceeds just as smoothly. Substitution of the energy (17) into 7c, . (10) 
yields 
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which result concurs with the direct solution of (14) f or the separate 
turning points 
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Corresponding agreement does not prevail, on the other band, between the 
f functions of a. rotating Morse isolccule. 


It is thus anticipated that a full-scale SKR treatment (f and » both) 
will not be so successful as (19) find (20) because cf the somewhat artificial 

cuacept •£ i vu tit a ia cruiiutii d through , Eq. (11); more precisely, since 
the WkB energy condition (i) can reveal at most the different «/ between 
classical turning points at a specified level, Klein’s g function must be 


regarded as a clever artifice for extending the theory and providing additional 

(but i»s# reliable) intonation. Failure of Eqs. (10) and (11) to attain 

complete accuracy for r and r . individual lv for a Morse potential 

sax mm 

cannot, therefore, be attributed to any failure per se of the sh? method, 
doe of our goals, then, is to eliminate as nonossential all reference to ro- 
tation in analyses of nonrotating diatowic ssolecules. 

AH ALTERNATE FORMOLATIOH 


As a substitute for the above introduction of rotational parameters , 
we now propose a compromise between the RjCR ideal of no predetermined an- 
alytic potential and those formulations that do so begin. Since sost 
analytic potentials, the Morse included despite its obvious deficiencies 


10 - 


at small separations, are quite adequate over considerable portions of 

their inner or repulsive legs, it seems reasonable to apply such a function 

to those regions only and subsequently ear-line with the SKR fonwila (10) for 

determination of the attractive part. Good repulsive branches also are avail 

able in many instances froa the reduction of high energy scattering data. 

Our second simplification is related, though rather loosely because of 

their retention of the g function, to an approximate technique tievised by 
IS 

Singh and Jain. A critical analysis of that aetbod is therefore desirable 
and we start with the diagram 


r *.x >B id y >l 
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-CD- 


-> f y (V) 


( 5 ) 


6* 

v 




■*» (RKR) 


(D 


(<:>) 


v 


m 


~ (5) > ^ CV„) . 


where the subscripts and superscripts a and m refer to analytic and Horse 
respectively. 

The various paths ?re explained as follows; 

Path (I) : The only argument stated explicitly by Singh and Jain re- 

duces to the fact that a good description of V, though perhaps only within 
a very narrow Av about the chosen quantum number, trivially guarantees good 
f values in that range. 
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Path (2) : Their actual procedure , on: the other jjami, is to start with 

goo o vibrational energies, as depicted by sop.e analytic G* within the 
3s»e v lirits, and then apply Klein’s f integral to compute (MR). 

Path (3): Unless, of course, values of 6 thnt are guaranteed accept- 

able only within Av automatically ensure ce .rrcspcnJiagly acceptable 
V[r (v)J , the argument for path (1) is entirely irrelevant and «,l- 

ditional proof for path (2) i s required, such evidence is discussed. 

Pstn (2) continued: It is desirable, therefore, to establish ?. con- 

Tineine reason f&t direct eotmlinz between accurate eaerzy levels laid reliable 
f (RKR) values, all wit in son* specified Av ^nc? such argument , 

though very' crude aah in apparent conflict with t’.o indication fxora Eq. (7) 
that lower levels cannot be ignored, develops cror the fact f-at a substantia 
p%tio& of Klein* s f intagrsl is contributed by levels fairly close to the 
chosen energy. Accordingly, that neighborhood is by far ch*. best narrow 
region in which to require a ,-ood description of f! . hinco Horse energies 
are seldom very poor, even for a single application over the es.tire range 
of v , snd sines the ; rcatest errors accompany in? s fit of fq. (17) to Av 
nearly always arise where they count the least, ve ray hopefully assert a 
reasonable accuracy for the ccrrespon d**:- evsinstinn of f at any level. 

Path (4): If the intention of Singh Jain Is to sirrnlify the f?or© 

rigorous calculation, path (2) canvot V»c apj* Hoa conveniently to :i least 
squares fit of G with a •'olvnowiul higher than second fopre-c. Path f 4) 

V 
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Rwrely associates this quadratic with the Morse potential (a * m) 
best do fines it according to Hqs. (14. (17), an-i (IS), wnereas the 


that 

equilibrium 


separation is provided fttss 


r 

e 
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( 21 ) 


Path {$) : The equation V (r . ) ■ C 1 * is solved directly to yield, 

v * * sr max 5 «ir. v ' 1 * 

as in the case of ( 20 ) above, the quantity f® (V ). Sinc« the latter sseci- 
£i** twsytKf .lit* thfee Marsa pwrasseters, only r is required to reverse this 

path back towerd V . 

1 s 

Path (6): Because of the property of first-oreer 4KB theory to repro- 

duce exactly a Morse potential, it is clear that r*“ (PKi.) 5 ^ (V^) ; con sc- - 
fwatl/, these is ns advantage to solving Klein’s f integral with 0£. 

What Singh and Jain evidently fail to recognize is that quadratic energy 
levels arc synonosvms with s Morse potential. Their sequence of path (2) 
plus the introduction of Klein’s g function is identical, exclusive of 


additional appro xiisat ions related to rotation, with the much simpler and far 
sore direct route of path (4), 4s an a least equivalent replaceaent for the 
Singh and Jain fumulation, we therefore adept the follewinp expressions: 

1/2 «»2 (mx<f) 1/2 i 
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and 


* w(v+|) - wx(v+~) 


(24) 


where w and «x sre determined from a least squares fit of the four ob- 
served vibrational levels closest to the one chosen. One need not, of course, 
both - with f at all, but it is a convenient quantity for comparison with 
rigorous MR calculations. 

CALCULATIONS ON fUfX 1 ! *) 

- 5 

Ground-state hydrogen is especially suited for a check of the afeove- 

aethod for three reasons: In the first place, Morse's potential is known to 

be a poorer fit for this than for Rost diatomic rtolecules; secondly, reliable 

|0 

data is available up to the dissociation 1 isit ; thirdly, nefssF.an, et al, 
have provided good RKR turning points. In order to v.ake the test even more 
stringent, we further single out the fourteenth vibration level (v * 11) , 
which energy is roughly 93% of that of dissociation and therefore lies 
in the region least adaptable to an overall analytic curve-fit and most 
subject to the errors discussed in the preceding section. 

Morse potentials were employee’ in several different ways in the present 


calculation. The third and fifth columns of Table I were obtained from known 


values of D and G , each potential extending over the entire range of 
v but with the saplicaticns differing according to whether G, _ or (T. „ 
(spectroscopic data) was used in tqs. (19) and (20), and the spectroscopic 
constants were computed from 


m » 4L< I 1-(1-D‘ 1 G ) 1/2 
v o J 


*x = 40 1 l-Cl-D^C ) x/m \ 

j ° J 


i 2 


(25) 


Colunn four also represents a single potential over the entire range, only 


this tisao with the spectroscopic constants determined froft n least squares 

fit of all 1& observed i avals. Again Gj S was used and D was dadncad 

frost Eq. (18) . 


As expected, the best results are associated with the last colusm, 
where the procedure was identical with the fourth except for the use only of 


the top four energies in the least squares determination of W and wx. 


Trough they differ slir-htly (within 1%) fro**- these of Tobias 
heissnaa's calculations were regarded as exact in nil cases. 


. . 1 
an.. sneers nee, 

it is interestin 


tc note that despite the fact that a better potential would have yielded s 
smaller error in r . , that discrepancy is in the right direction to 
actually increase the reliability of r ^ from Eq. (2)1 the same world 
have been accomplished automatically, however, fro» the r form of (20) 

Hi &X 

, jj) 

xn conjunction with the sane C. . and independently of r , . 


SUMMARY 


Two new not hods. one of which involves a transformation that simul- 
taneously renoves the upper -licit singularities of Klein integrands and 
establishes a single integration variable through out, have been introduced 
for finding eiatoaic potentials more or less directly from spectroscopic 
data. The ether is an alternate procedure that effectively replaces the 
Ml integrals* while at the s*m tine deriving a sort of semi -rigorous 
support ftws that foriulation. Janaerical calculations on normal hydrogen 
revealed s discrepancy of less than Q.S % for the attractive branch of the 
potential; accordingly, the alternate technique, together with Keitier- 
London type approximations at large separations,"'' should provide an 
adequate representation for nany collision problem* - including those 
conerraed with transport ana other nonequilibrium phenomena in dilute gases, 
because of the ease of calculation without -electronic cctsmiters, the alternate 
tnstmnt sipht aptly be termed & -peer man’s KK” . ' 






Table 1. Classical turning points for 
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(v » --34) 
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0 , CM 

33, 274 C 

35,274 

40,025 

55,274 

38,895 

0 

i A A 

1.131 

1.116 

1.051 

1.358 

1.139 

* error 

A 

-1.33 

-7.07 

20.1 

0.71 

$c 

r . ,A 
mm 

0.413 

0.393 
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0.389 

0.408 

" error 


-4.84 
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-1.21 

o 

r , A 
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? 

4* « V i mf 

2.625 

2.48S 

3.104 

2.686 

9 error 

«> 

-1.87 

-6.99 

16.0 

0.41 
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0 , where D - D ♦ C , was taker: from G, Herzbere, Snectra of diatomic 

e* co 1 

Molecules (D. Van Hostrrand, Inc., Princeton, h'ev Jersey, 1950), 2nd ad., 
p. 530. 
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